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The real-time dynamics of photoexcited and photoionized rubidium (Rb) atoms attached to helium 
(He) nanodroplets is studied by femtosecond pump-probe mass spectrometry. While excited Rb 
atoms in the perturbed 6p-state (Rb*) desorb off the He droplets, Rb^ photoions tend to sink into 
the droplet interior when created near the droplet surface. The transition from Rb^ solvation to full 
Rb* desorption is found to occur at a delay time r ~ 600 fs for Rb* in the 6pE-state and r ~ 1200 
fs for the GpH-state. Rb^He ions are found to be created by directly exciting bound Rb*He exciplex 
states as well as by populating bound Rb^He-states in an photoassociative ionization process. 


I. INTRODUCTION 

The laser-induced dynamics of pure and doped helium 
(He) nanodroplets is currently attracting considerable at¬ 
tention [1-7]. While the superfluidity of He nanodroplets 
has been tested in numerous key experiments by prob¬ 
ing stationary properties [8-10], the impact of the quan¬ 
tum nature of the droplets on their dynamic response 
to impulsive excitation or ionization is much less well 
established. As a prominent recent example, the rota¬ 
tional dynamics of various molecules embedded in He 
droplets induced by impulsive alignment was found to 
be significantly slowed down and rotational recurrences 
were completely absent [5]. This indicates that substan¬ 
tial transient system-bath interactions are present during 
the laser pulse. In contrast, the vibrational dynamics of 
rubidium (Rb) molecules Rb 2 attached to the surface of 
He nanodroplets revealed only slow relaxation and de¬ 
phasing proceeding on a nanosecond time scale [2, 11]. 

Various recent experimental and theoretical studies 
have addressed the dynamics of solvation and desolva¬ 
tion of ionized or excited metal atoms off the surface of 
He nanodroplets [7, 12-21]. So far, these studies have 
concentrated on measuring the total yield and the final 
velocity of the ejected atoms as a function of the atomic 
species and the electronic state of excitation. In this pa¬ 
per we present the first time-resolved characterization of 
the desorption process of Rb atoms off the surface of He 
nanodroplets upon excitation to the droplet-perturbed 
states correlating to the 6p atomic orbital. The experi¬ 
mental scheme we apply is femtosecond (fs) pump-probe 
photoionization in combination with time-of-fiight mass- 
spectrometry. We find that the yield of detected Rb+ 
photoions as a function of delay time r between the ex¬ 
citing pump and the ionizing probe pulses is determined 
by the interplay of the repulsive interaction of excited 
Rb* with respect to the He surface and the attractive 
interaction of the Rb+ ion with the He surface induced 
by photoionization. 

The Rb*-He droplet repulsion initiates the desorption 
of the Rb* atom off the He droplet surface. Except for the 
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lowest excited state of Rb, 5pi/2, all excited states up to 
high Rydberg levels experience strong repulsion from He 
droplets [22, 23]. In contrast, the Rb+-He droplet attrac¬ 
tion causes the Rb+ ion to fall back into the He droplet 
when created near the He droplet surface at short delay 
times [19, 21]. Atomic cations are known to form stable 
“snowball” structures consisting of a cationic core which 
is surrounded by a high density shell of He atoms. As a 
result, free Rb+ ions appear in the mass spectrum only 
after a characteristic pump-probe delay time which 
depends on the state the Rb atom is initially excited to. 

In addition to neat Rb+ atomic ions, the photoioniza¬ 
tion mass spectra contain Rb+He molecular ions in the 
full range of laser wavelengths correlating to the droplet- 
perturbed Rb 6p-state. The occurrence of such molecu¬ 
lar ions has previously been interpreted by the formation 
of metastable ‘exciplex’ molecules [1, 7, 14, 16, 24, 25]. 
These bound states of excited metal atoms and one or 
few He atoms can be populated either by a tunneling 
process [12, 26-28] or by direct laser-excitation of bound 
states in the metal atom-He pair potential [14, 16, 17, 29]. 
In the former case, exciplex formation times > 50 ps are 
expected [24, 27], whereas in the latter case, exciplexes 
are created instantaneously. Thus, previous pump-probe 
measurements revealing exciplex formation times of 8.5 
and 11.6 ps for Rb^He and Rb^He, respectively, upon ex¬ 
citation into the droplet-perturbed 5p3/2-state could not 
be consistently interpreted [24]. 

In the present study we observe a time-delayed increase 
of the Rb+He signal as for Rb+ indicating that the pump- 
probe dynamics is primarily determined by the competi¬ 
tion between desorption of the Rb*He exciplex off the He 
droplet surface and the Rb+He cation falling back into 
the He droplet interior. Moreover, a pronounced max¬ 
imum in the Rb+He signal transients indicates that an 
additional Rb+He formation channel besides photoion¬ 
ization of Rb+He exciplexes is active - photoassociative 
ionization (PAI) of the desorbing Rb atom and a He atom 
out of the droplet surface. PAI is a well-known process 
where a bound cationic molecule or complex is formed 
by photoionization or photoexcitation into autoionizing 
states of an atom or molecule of a collision complex [30]. 
PAI is a special case of traditional associative ionization 
where a bound molecular cation is formed in a binary 
collision of an electronically excited atom [31]. In either 
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case the binding energy is taken away by the electron 
emitted in the process. 


II. EXPERIMENTAL SETUP 

The experimental setup is similar to the previously 
used arrangement [11, 14] except for the ionization and 
detection schemes. He droplets are produced by a con¬ 
tinuous supersonic expansion of He 6.0 through a 5 /im 
nozzle at a pressure of 50 bar. The transversal velocity 
spread of the beam is reduced by placing a 400 fim skim¬ 
mer 13 mm behind the nozzle. Unless otherwise stated, 
the nozzle temperature is kept at 17 K. This results in 
a log-normal distribution of the He droplet size with a 
mean size of 1.1 x 10^ He atoms. Subsequently, the 
droplet beam passes a mechanical chopper and a Rb-filled 
cell of length 1 cm, stabilized at a temperature of 85 °C. 
At the corresponding vapor pressure, most droplets pick 
up on average one Rb atom following poissonian statis¬ 
tics. By overlapping the droplet beam with the output of 
the fs laser, we resonantly excite and ionize the dopant 
atom. 

In contrast to previous studies, we use amplified fs laser 
pulses generated by a regenerative amplifier operated at 
a pulse repetition rate of 5 kHz. At this repetition rate, 
multiple excitations of Rb atoms by subsequent pulses 
from the pulse train are safely excluded. The pulses are 
frequency-doubled in a BBO crystal resulting in a pulse 
duration of tp = 120 fs with a variation for different laser 
center wavelengths of 20 fs. Two identical, time-delayed 
pump and probe pulses are generated by means of a me¬ 
chanical delay line. The laser beam is focused into the 
vacuum chamber using a 30 cm lens which leads to a peak 
intensity in the range of 5 x 10^^ Wcm“^. 

Photoions are detected by a time-of-flight (TOP) mass 
spectrometer in Wiley-McLaren configuration mounted 
in-line with the He droplet beam [32]. At the end of the 
drift tube a high negative potential is applied to further 
accelerate the arriving ions which boosts the efficiency 
of detecting large cluster masses in the 10^ amu range 
using a Daly-type detector [33]. The latter consists of 
a Faraday cup, a scintillator with an optical bandpass 
interference filter and a photomultiplier tube. In case 
of electron detection, a simple electrode setup consisting 
of a repeller, an extractor grid and a channeltron detec¬ 
tor with positive entrance potential is used. For both 
detectors, the resulting pulses are amplified, threshold- 
discriminated and acquired by a fast digitizer. When 
detecting heavy masses a counting unit is used. 


III. Rb DESORPTION DYNAMICS 

In the present paper we concentrate on the fs pump- 
probe dynamics of Rb atoms attached to He nanodroplets 
which are excited to droplet-perturbed states correlating 
to the atomic 6p-state. These states have previously been 



FIG. 1. Potential energy diagram of the Rb-He nanodroplet 
complex. Vertical arrows depict the photo-excitation and ion¬ 
ization processes. The potential curves of the neutral Rb- 
He 2 ooo complex are taken from [23], the one of the Rb^-He 2 ooo 
complex is obtained from the Rb^-He pair potential [36] on 
the basis of the He density distribution of the groundstate 
RbHe 2 ooo complex [37]. The two peaks plotted vertically on 
the left-hand scale show the expected excitation spectrum 
based on these potentials. 

studied using nanosecond pulsed excitation and velocity- 
map imaging of photoions and electrons [14, 16]. Due 
to the interaction of the excited Rb atom with the He 
droplet surface, the 6p-state splits up into the two states 
6pE and bpH according to the pseudo-diatomic model 
which treats the whole He droplet. He at, as one con¬ 
stituent atom of the RbHeAr complex [23, 34, 35]. 

Using the RbHeAr pseudo-diatomic potential curves for 
the 5sE electronic groundstate and the 6pE, H excited 
states we compute the Franck-Condon profiles for the 
expected vertical excitation probability using R. LeRoy’s 
program BCONT [38]. The corresponding transition 
probability profile is depicted on the left-hand side of 
Fig. 1. The experimental excitation spectrum is in good 
agreement with the calculated one apart from the fact 
that the experimental peaks are somewhat broader [14]. 
Since both 6pE and bpH pseudo-diatomic potentials are 
shifted up in energy by up to 1200 cm“^ with respect to 
the atomic 6p level energy, we expect strong repulsion 
and therefore fast desorption of the Rb atom off the He 
droplet surface to occur following the excitation. 

However, upon ionization of the excited Rb atom by 
absorption of a second photon (vertical arrow on the 
right-hand side of Fig. 1), the interaction potential sud- 
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FIG. 3. (a) Typical mass spectra recorded for Rb-doped He 
nanodroplets by fs photoionization taken at a center wave¬ 
length A = 415 nm and a 5 ps pump probe delay. In addition 
to the atomic isotopes ®^Rb^ and ®^Rb^ the mass spectra 
contain Rb^He and Rb^He 2 molecular ions, (b) An extended 
view of mass spectra taken at various nozzle temperatures us¬ 
ing single fs pulses at A = 415 nm reveals the presence of large 
masses of unfragmented ion-doped He droplets Rb^Hew- 


FIG. 2. Pump-probe transient Rb^ ion count rates recorded 
for various wavelengths A of the fs laser. At A > 409 nm, 
excitation occurs predominantly to the GpH-state, at A < 409 
nm predominantly the 6pE-state is excited. The thin smooth 
lines are fits to the data. 


denly turns weakly attractive. Thus, the Rb+ ion may 
be expected to turn around and to fall back into the He 
droplet provided ionization occurs at short delay times 
after excitation such that the desorbing Rb* picks up 
only little kinetic energy 


Ekin,Rh*{R) < Epot,nh+{R)- ( 1 ) 


Here, Rb+ (^) denotes the lowering of the potential 
energy of the Rb+ ion due to the attractive interaction 
with the He droplet at the distance R from the droplet 
surface. Eq. 1 holds for short distances R < Rc falling 
below a critical value Rc. When assuming classical mo¬ 
tion, we can infer from Eq. 1 the critical distance Rc for 
the turn-over. Erom simulating the classical trajectory 
R{t) we can then obtain the delay time Tc at which the 
turn-over occurs. In the following we refer to Tc as ‘fall¬ 
back time’. Thus, when measuring the number of free 
Rb+ ions emitted from the He droplets by pump-probe 
photoionization we may expect vanishing count rates at 
short delays r < Tc due to the Rb+ ions falling back into 
the droplets, followed by a steep increase and subsequent 
constant level of the Rb+ signal at delays r > Tc. 


A. Experimental results 

Eig. 2 shows the transient Rb+ ion signals measured 
by integrating over the ^^Rb and ^^Rb mass peaks in 
the time-of-flight mass spectra recorded for each value 
of the pump-probe delay. The shown data are obtained 
by subtracting from the measured ion signals the sum of 
ion counts for pump and probe laser pulses only. The 
error bars stem from error propagation taking into ac¬ 
count the uncertainties associated with the different sig¬ 
nal contributions. By tuning the wavelength of the fs 
laser A we can excite predominantly the GpH (A > 409 
nm) or the 6pE-states (A < 409) of the ROHcat complex. 
As expected, we observe a step-like increase of the Rb+- 
yield at delays ranging from 600 fs (A = 401 nm) up to 
about 1500 fs (A = 415 nm). The signal increase occurs 
at shorter delays when exciting into the more repulsive 
6pE-state because the Rb atom moves away from the 
He droplet surface faster than when it is excited into the 
shallower GpH-state. The rising edge of the signal jump is 
extended over a delay period of about 400 fs, partly due 
to the finite length and bandwidth of the laser pulses. 
Desorption along the GpH-potential appears as an even 
smoother signal rise, indicating that a purely classical 
model is not suitable for reproducing the observed dy¬ 
namics. Eor laser wavelengths A < 409 nm we observe a 
weakly pronounced double-hump structure with maxima 
around 800 and 1800 fs, respectively, which we discuss in 
section IIIB. 

Before discussing our model calculations for these tran¬ 
sients, let us first examine the measured time-of-fiight 
mass spectra in more detail. Eig. 3 (a) depicts a represen¬ 
tative mass spectrum in the mass range around 100 amu 
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at a pump-probe delay of 5 ps and a center wavelength 
A = 415 nm. The spectrum is averaged over 5000 laser 
shots. Clearly, the dominant fragments in this mass 
range are neat Rb+ ions at 85 and 87 amu, where the 
different peak heights reflect the natural abundances of 
isotopes (72 and 28 %, respectively). Even when ion¬ 
izing with single laser pulses the mass spectra contain 
bare Rb+ ions at a low level. We attribute this to a 
fraction of the Rb atoms desorbing off the droplets and 
subsequently ionizing within the laser pulse. A contri¬ 
bution to the Rb+ signal may come from free Rb atoms 
accompanying the droplet beam as a consequence of the 
detachment of the Rb atom from the droplet during the 
pick-up process. Aside from neat Rb+ atomic ions, the 
pump-probe mass spectra feature peaks at 89, 91, and 
95 amu, which evidence the formation of Rb+He and 
Rb+He 2 molecular ions. These masses are usually at¬ 
tributed to photoionization of bound metastable Rb*He 
exciplexes [1, 14, 17, 24]. 

In addition to these discrete mass peaks, we mea¬ 
sure extended mass distributions reaching up to 64,000 
amu using our time-of-flight mass spectrometer which 
is optimized to detecting cluster ions. These distribu¬ 
tions are in good agreement with the size distributions 
of pure He nanodroplets generated in a sub-critical ex¬ 
pansion [39, 40]. From comparing the peak areas of the 
light masses Rb+, Rb+He^, n = 1,2 with those of the 
heavy droplet ions Rb+He 7 v we deduce that by ioniz¬ 
ing with single pulses a fraction of < 10% of the doped 
He droplets fragments into free atomic or molecular ions. 
The larger part of the ionized Rb-doped He droplets gen¬ 
erates unfragmented Rb+HcTv due to the sinking of the 
Rb+ ion into the He droplet and the formation of a sta¬ 
ble snowball complex [19]. When adding an additional 
time-delayed probe pulse we may expect to alter this ra¬ 
tio by depleting the unfragmented Rb+HcTv fraction in 
favor of creating free ions Rb+ and Rb+Hei ^2 ions after 
desorption. 

Indeed, the delay-dependent peak integrals of the mea¬ 
sured mass peaks at A = 400 nm confirm this picture, 
see Fig. 4 (a) and (b). While the atomic Rb+ ion sig¬ 
nal sharply increases around 600 fs and remains largely 
constant for longer delays, the Rb+Hciv signal displays 
the opposite behavior. The maximum signal level at zero 
delay significantly drops around r = 600 fs and remains 
low for long delay times. 

In addition to the mass-resolved ion signals we have 
measured the total yield of photoelectrons, depicted in 
Fig. 4 (c). From comparing the electron counts with 
and without blocking the He droplet beam we find that 
for pump-probe ionization >79% of photoelectrons cor¬ 
relates with the Rb-doped He droplet beam, < 21% is 
attributed to ionization of Rb and other species in the 
background gas. The observation that the electron count 
rate remains constant within the experimental scatter in 
the entire range of pump-probe delays indicates that the 
photoionization efficiency (cross-section) of a Rb atom is 
largely independent of its position with respect to the 
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FIG. 4. Transient ion and electron signals measured at A = 
400 nm. The ion signal traces (a and b) are obtained from 
integrating over the free atomic Rb^ ion peaks and over the 
charged He droplet mass distribution, respectively. The total 
photoelectron signal (c) is measured using a simple electron 
detector. The thin smooth lines are fits to the data. 

He droplet surface. These observations further support 
our interpretation of the step-like increase of Rb+ counts 
in terms of the competition between desorption of ex¬ 
cited Rb atoms and solvation of Rb+ cations into the He 
droplets. 

Fig. 5 (b) displays a compilation of the critical delays 
for all measured laser wavelengths which we obtain by 
fitting the experimental data with an error function, 

fRb+ (i) = ^ ■ {erf [{t - rc)/(T] 1} (2) 

of variable amplitude A, width a and position Tc. Shown 
are the results for the raw Rb+ and Rb+Hew transients 
as well as those obtained by fitting the sum of the tran¬ 
sients of Rb+ atomic and Rb+He molecular ions. In par¬ 
ticular for the bpH state, the latter signal more accurately 
reflects the dynamics of the fail-back process than the 
individual Rb+ and Rb+He transients since additional 
transient redistribution of population between Rb+ and 
Rb+He channels, which we discuss below, cancels out. 
Correspondingly, the fitted time constants of the summed 
Rb+ and Rb+He transients and those of Rb+Hcw are in 
good agreement. This confirms our conception that the 
light ions fall back to produce heavy cluster ions at short 
delays. Fig. 5 (c) will be discussed in section IV. 


B. Simulations 

Further support for our interpretation of the experi¬ 
mental findings is provided by classical trajectory simula¬ 
tions of the dynamics of the pump-probe process. In this 
model, the Rb atom and the He droplet surface are taken 
as two point-like particles which propagate classically ac¬ 
cording to the pseudo-diatomic model potentials [23, 37]. 
Note that these potentials were calculated based on the 
minimum-energy configuration of a droplet consisting of 
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FIG. 5. Simulated (a) and experimental (b) fall-back times 
as a function of the laser wavelength, derived from the ris¬ 
ing edges of the pump-probe transients. The curves in (a) 
are obtained for various effective masses mnen — ^eff^He of 
the He droplet in units of the He atomic mass mne- The dif¬ 
ferent symbols in (b) denote the experimental fit results for 
Rb^, Rb^He+Rb^ and Rb^Hew signals. Panel (c) shows 
the exponential decay constants from fits of the Rb^He ion 
transients with Eq. (4). 


N = 2000 He atoms subjected to the external potential of 
an Rb atom in the electronic ground state. The classical 
equation of motion 




dV{R) 

dR 


(3) 


is solved numerically. Here, V = jj Rb+(^) denotes 
the potential curves of the excited and ionic states, and 
R{t) is the distance between the Rb atom and the He dim¬ 
ple at the droplet surface. The initial value of the Rb-He 
droplet distance is the position of the minimum of the 
groundstate potential well (6.4 A). Eq. 3 is first solved 
for the neutral excited state potential or Vu up to the 
pump-probe delay time r. Subsequently, the Rb atom is 
considered to be ionized and the particle is propagated 
further using the ionic Rb+-HeAr potential The 

reduced mass /i = mRbUiHen/(^Rb + is given by 

the mass of the Rb atom or ion, ruRb, and the effective 
mass of the He droplet, We set mRen = 40 amu 

for the propagation of the excited as well as for the sub¬ 
sequent propagation of the Rb+ ion with respect to the 
He droplet. This value is based on previous experimental 
as well as theoretical findings [16]. 

The motion of the excited and subsequently ionized 
Rb atom with respect to the He droplet surface is illus¬ 
trated in Fig. 6 for different initial conditions. The time- 



FIG. 6. Glassical trajectories of the excited and ionized Rb 
atom initially located in a dimple near the He droplet surface, 
(a) At bpH-excitation (A = 415 nm) and long pump-probe 
delay r = 500 fs the Rb atom fully desorbs off the He droplet 
and propagates as a free Rb^ cation after ionization, (b) At 
shorter r = 400 fs the Rb^ ion turns over and falls back 
into the He droplet. The schemes at the bottom visualize the 
dynamics for increasing time from left to right. 


dependent positions of the Rb atom and the He surface 
are depicted as red and blue lines in the upper parts. The 
lower parts are graphical visualizations of the dynamics. 
Fig. 6 (a) depicts the case when the excitation of the Rb 
atom, which is initially located in the groundstate equi¬ 
librium configuration of the RbHcAr complex, occurs at 
t = 0 and ionization is delayed to r = 500 fs. The laser 
wavelength is set to A = 415 nm where the motion fol¬ 
lows the OpH-potential. In this case the excited Rb atom 
fully desorbs off the He droplet and continues to move 
away from the droplet after its conversion into an ion. In 
the case of shorter delay r = 400 fs between excitation 
and ionization, shown in Fig. 6 (b), the Rb atom turns 
over upon ionization as a result of Rb+-Heiv attraction 
and falls back into the He droplet. 

For assessing the effect of an initial spread of Rb-He at 
droplet distances R due to the broad laser bandwidth 
and of the finite length of the laser pulses tp we extend 
the classical trajectory calculation to a mixed quantum- 
classical simulation which includes an approximate de¬ 
scription of the quantum wave packet dynamics of the 
system. 

The initial wave packet is obtained by transforming the 
spectral profile of the laser into a distribution as a func¬ 
tion of R using the potential energy difference between 
the initial 5sE and the final 6pE, H pseudo-diatomic 
states. We use a Gaussian-shaped laser profile with a 
full width at half maximum, Az/, inferred from measured 
spectra. Typically Au ^ 2 nm, depending on the center 
wavelength of the laser. This corresponds to the instan¬ 
taneous creation of a wave packet in the excited state 
centered around the minimum of the groundstate poten¬ 
tial. 

For simulating the dynamics the wave packet is approx¬ 
imated by 25 segments i and each segment is propagated 
individually according to Eq. 3 where R{t) is replaced 
by Ri{t) representing the Rb-He at distance for the i-th 
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FIG. 7. Semiclassical simulations of the yield of free Rb^ ions 
created by excitation and time-delayed ionization for various 
center wavelengths A of the laser pulses. See text for details. 


segment. Convergence of the final results with respect to 
the number of segments has been checked. 

This simplified description of the wave packet dynam¬ 
ics is justified because no quantum interference effects are 
expected for this simple dissociation reaction. Compar¬ 
ison with the full quantum simulation of the desorption 
process yields excellent agreement within the propaga¬ 
tion time range relevant to the experiment. 

Simulated transient yield curves as a function of the 
pump-probe delay r are obtained by taking the weighted 
sum of the segments which have propagated outwards 
up to very large distances after long propagation times. 
This sum we identify with the fraction of desorbed atoms. 
Those segments which have turned over towards short 
distances are considered to contribute to the Rb+ ions 
falling back into the droplet. For those segments the 
condition formulated initially (inequality (1)) is fulfilled 
implicitly. The finite duration of the excitation process 
is taken into account by convolving the resulting yield 
curves with the autocorrelation function of the two laser 
pulses. 

The resulting simulated yields of free Rb+ ions as a 
function of r are depicted in Fig. 7 for various center 
wavelengths A of the laser pulses. The obtained curves 
qualitatively resemble the experimental ones in many re¬ 
spects. Excitation at long wavelengths A > 409 nm, 
at which predominantly the more weakly repulsive bpll 
state is populated, induces a smooth signal increase at 
about T = 400 fs. At A < 409 nm, where predomi¬ 


nantly the 6pE-state is excited, the signal rise occurs 
around r = 210 fs, considerably earlier than for bpll 
excitation. This result qualitatively agrees with the ex¬ 
perimental finding see Fig. 5 (b). Moreover, the superpo¬ 
sition of the two rising edges at intermediate wavelengths 
A ^ 409 nm may provide an explanation for the double¬ 
hump structure observed in the experimental Rb+ tran¬ 
sients at A < 409 nm. However, the simulated rising 
edges occur at significantly shorter delay times than in 
the experiment, roughly by a factor 2 for excitations to 
the 6pE-state and up to a factor 4 for the bpH-state. 

The discrepancy between the experimental results and 
those of the simulations, shown in Fig. 5 (b), is present 
even when assuming very large effective masses of the 
He droplet > 1000 amu. We attribute this dis¬ 

crepancy to the limited validity of our model assump¬ 
tions. In particular the interaction potentials we use 
were obtained on the basis of the frozen He density dis¬ 
tribution for the RbHcAT groundstate equilibrium con¬ 
figuration [23, 37]. However, transient deformations of 
the He droplet surface in the course of the dynamics 
are likely to significantly modify the effective Rb*-HeAr 
interactions. Recent time-dependent density functional 
simulations show a complex ultrafast response of the 
He droplet to the presence of a Rb+ ion near the sur¬ 
face [21]. In particular when the desorption dynamics 
is slow (H-state) a complex reorganization of the He 
droplet surface during the Rb desorption process may 
be expected [16]. A clear manifestation of the break¬ 
down of the simple pseudo-diatomic model is the forma¬ 
tion of Rb*He exciplexes which we discuss in the follow¬ 
ing section. Recently, M. Drabbels and coworkers sug¬ 
gested that the pseudo-diatomic potentials of the excited 
Na*HeAr complex may be transiently shifted and even in¬ 
tersect [17, 28]. Detailed three-dimensional simulations 
including the full spectrum of properties of He droplets 
are needed to provide an accurate description of this kind 
of dynamics [16, 20, 21, 41]. Experimentally, the time 
evolution of the interaction potential energies will be vi¬ 
sualized by means of fs time-resolved photoelectron spec¬ 
troscopy in the near future. 


IV. RbHe+ DYNAMICS 

Aside from free Rb+ ions, fs photoionization of Rb- 
doped He nanodroplets generates Rb+He^, n = 1,2 
molecular ions. Relative abundances reach up to 31% 
and 1.5%, respectively, measured at A = 415 nm corre¬ 
sponding to the 6pn excitation, see Eig. 3. At A = 399 
nm (6pE-excitation), abundances are 4% and 1%, respec¬ 
tively. Eree Rb+He ions are associated with bound states 
in the Rb*He excited states pair potentials, so called ex¬ 
ciplexes. Both the 6pE and the bpH-states of the RbHe 
diatom feature potential wells which sustain bound vi¬ 
brational states that can be directly populated by laser 
excitation out of the groundstate of the RbHe at com¬ 
plex [14, 29]. Thus, exciplexes are directly created in 

















7 


1.0 


0.5 


0.0 

1.0 


0.5 


3 0.0 
^ 1.0 

0.5 

c 

g) 

0.0 
X 1-0 

cc 

0.5 


0.0 

1.0 

0.5 


0.0 

0 1000 2000 3000 4000 5000 

Pump-probe delay (fs) 

FIG. 8. Experimental yields of Rb^He molecular ions as a 
function of pump-probe delay for various center wavelengths 
of the laser pulses. The thin smooth lines are hts to the data. 
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a process akin to photoassociation, in contrast to previ¬ 
ously observed Na*He and K*He exciplexes which were 
formed by an indirect tunneling process upon excitation 
of the lowest pll-states [27, 28]. 

Exciplex formation is the only route to producing 
Rb+He ions by photoionization using continuous-wave or 
nanosecond lasers, where ionization takes place at long 
delay times when the dynamics of exciplex formation and 
desorption off the droplets is long complete. In fs ex¬ 
periments, however, ionization can be triggered before 
or during the process of desorption of the excited atom 
or exciplex off the droplet surface. In this case, due to 
the attractive Rb+-He potential a bound Rb+He molec¬ 
ular ion can be formed upon ionization, even if the ex¬ 
cited Rb*-He interaction does not sustain bound states 
of the neutral diatom. The process of inducing a molec¬ 
ular bond between two initially unbound neutral species 
by photoionization is known as photoassociative ioniza¬ 
tion [31]. 


A. Experimental results 

The transient yield of Rb+He for various laser wave¬ 
lengths is displayed in Fig. 8. Similarly to the Rb+ tran¬ 
sients, we measure vanishing Rb+He pump-probe signal 
contrast around zero delay. For increasing laser wave¬ 
length from A = 399 up to 418 nm, which corresponds to 



FIG. 9. Glassical trajectories of the Rb-He-Hciv three-body 
system at A = 409 nm. The schemes at the bottom visualize 
the various dynamics for increasing time from left to right, (a) 
The excited Rb atom departing from the He droplet surface 
suddenly experiences Rb^He pair attraction upon ionization 
at T = 400 fs (a). Gonsequently, a He atom attaches to the Rb 
atom while it leaves the droplet, (b) For short delay r = 350 fs 
at A = 409 nm a Rb^He molecule forms as in (a) but the 
attraction towards the He droplet makes it turn over and fall 
back. 


the crossover from the 6pE to the GpH excited pseudo- 
diatomic states, a step-like increase of the Rb+He ion 
signal occurs at delays ranging from r = 500 fs up to 
about 2000 fs. Besides, at A < 415 nm we measure a 
transient overshoot of the Rb+He signal by up to about 
100% of the signal level at long delays. The transient 
yield of Rb+He is fitted using the model function 

fRb+He{t) = fRb+{t){Ee -h 1). (4) 

As for the Rb+ case, /^ 5 + (t) models the fall-back dynam¬ 
ics by Eq. 2. Additionally, the exponential function with 
amplitude E and time constant te takes the transient 
overshoot into account, whereas the additive constant 
account for a r-independent Rb+He formation channel. 
The exponential time constants te are plotted as black 
circles in Fig. 5 (c). To obtain these values, the param¬ 
eters Tc and a are taken as constants from the fit of the 
sum of Rb+ and Rb+He signals with Eq. 2. Here we 
make the assumption that the fail-back dynamics is only 
weakly perturbed by the attachment of a He atom to the 
Rb atom or ion, which is confirmed by our simulations. 


B. Simulation of photoassociative ionization 

For a more quantitative interpretation of the Rb+He 
transients we extend our classical and mixed quantum- 
classical models to the one-dimensional three-body 
problem in the second stage of the calculation after ion¬ 
ization has occurred by including one individual He atom 
out of the surface layer. The classical trajectories are now 
obtained by solving three individual coupled equations 
of motion for the three individual particles Rb+, He, and 
He^. The Rb*-He 7 v interaction leading to desorption is 
represented by the pseudodiatomic potentials as before. 

























overlap integrals 


2 



Pump-probe delay (fs) 

FIG. 10. Simulations of the yield of free Rb^He-molecules 
created by excitation and time-delayed ionization for various 
center wavelengths of the laser pulses. See text for details. 


The Rb+-He dynamics is initialized by the velocity and 
distance of the dissociating Rb*HeAr complex at the 
moment of ionization. The Rb+-He pair interaction is 
given by the Rb+-He pair potential [36] augmented by a 
16.7 cm“^ deep potential step to account for the He-He at 
extraction energy as suggested by Reho et al. [14, 24, 27]. 

Exemplary trajectories are shown in Fig. 9 for two 
cases at A = 409 nm. For long pump-probe delays 
the Rb+ ion leaves the He droplet without attaching 
a He-atom, as shown in Fig. 6. However, there is a 
range of delays in which the desorbing Rb atom is far 
enough away from the droplet so that it will not fall 
back upon ionization, but it is still close enough to 
attract a He atom out of the droplet surface so as to 
form a bound molecular ion by PAI (Fig. 9 (a)). Fig. 9 
(b) illustrates the dynamics at short delay when the 
attractive forces acting between the Rb+ ion and the 
droplet surface prevent the full desorption and Rb+-He 
pairwise attraction leads to the formation of Rb+He. 

For simulating the transient Rb+He yields to compare 
with the experimental data shown in Fig. 8 we extend 
the mixed quantum-classical model for the desorption dy¬ 
namics of bare Rb described in Sec. HI B. It is augmented 
by computing the probability of populating bound vibra¬ 
tional states of the Rb+He molecule for each segment of 
the Rb wave packet upon ionization as the sum of spatial 


pOO 

= F / fpi{R, t) dR 

^ oo 


(5) 


Here t/;* denotes the i-th wave packet segment and (py 
stands for the vibrational wave functions of Rb+He cal¬ 
culated using R. J. LeRoy’s LEVEL program [42] for the 
6pF, H pair potentials of Rb+He [29]. The identifica¬ 
tion of bound and free Rb+He ions in the simulation is 
based on analyzing the final Rb+-He and Rb+He-HcAr 
distances after long delays r > 10 ps, respectively. The 
final probability P of detecting a Rb+He molecule is ob¬ 
tained by summing up the detection probabilities for ev¬ 
ery segment, 

P{r) = J2P?ir)-{pr^\r)+Pe.). ( 6 ) 

i 

In agreement with Eq. 4, denotes the desorption 
probability and Pex is the probability of creating a bound 
neutral Rb*He exciplex, which is assumed to occur in¬ 
stantaneously upon laser excitation and thus does not 
depend on r. Since the relative contributions of PAI 
and direct exciplex formation are not precisely known we 
resign from quantitatively modeling the relative efficien¬ 
cies of the two pathways leading to free Rb+He. Instead 
we adjust them to the experimental transients by taking 
Pex as a free fit parameter. The transient signal P{r) 
is finally convoluted with the intensity autocorrelation 
function of the laser pulses, as for the Rb+ transients. 

The resulting simulated yields of free Rb+He molecu¬ 
lar ions are depicted in Fig. 10. Clearly, the same gen¬ 
eral trends as for neat Rb+ ions are recovered: (i) at 
short delay times r < 200 fs the appearance of Rb+He 
is suppressed due to the falling back of the ion into the 
He droplet; (ii) longer laser wavelengths A > 409 nm 
(OpH-excitation) lead to weaker repulsion and therefore 
to the delayed appearance of free ions as compared to 
6pF-excitation at A < 409 nm. These results again qual¬ 
itatively agree with the experimental findings but the 
simulated appearance times are shorter by a factor 2-4, 
as shown in Fig. 5. As in the Rb+ case we attribute 
these deviations to the use of pseudo-diatomic potentials 
calculated for the frozen RbHcAr groundstate complex. 

Moreover, the simulation reproduces a signal overshoot 
around r = 300 fs at short wavelengths, which is due 
to the contribution of the photoassociative ionization 
channel. Association of a bound Rb+He ion is possi¬ 
ble only at sufficiently short Rb-He distances given at 
delays r < 600 fs for 6pF-excitation and r < 900 fs for 
OpH-excitation, respectively. At these delay times, the 
PAI signal adds to the signal due to ionization of Rb*He 
exciplexes formed directly by the pump pulse. Note that 
we have adjusted to the experimental curves the relative 
contributions to the total Rb+He signal arising from PAI 
and from exciplex ionization. Therefore our model does 
not permit a quantitative comparison with the experi¬ 
mentally measured signal amplitudes. A more detailed 
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three-dimensional simulation of the dynamics is needed 
for a fully quantitative interpretation [20]. Nevertheless, 
we take the simulation result as a clear indication that 
PAI is an additional channel producing He-containing 
ionic complexes which needs to be considered in exper¬ 
iments involving photoionization of dopants attached to 
He droplets when using ultrashort pulses. 

We note that in the particular case of exciting into the 
6pS-state of the RbHcTv complex, unusual Rb+He sig¬ 
nal transients may arise from the peculiar shape of the 
RbHe pair potential curve which features a local poten¬ 
tial maximum at intermediate Rb-He distance [14, 29]. 
This potential barrier causes the highest-lying Rb*He vi¬ 
brational states to be predissociative. Semi-classical es¬ 
timates yield predissociation time constants for the two 
highest vibrational levels v = b and 6 of 3.5 ns and 2.2 ps, 
respectively. However, these values significantly exceed 
the exponential decay times inferred from the measured 
transients (see Fig. 5 (c)). Moreover, we may expect 
that not only the highest vibrational levels are populated. 
Note that for the case of Na*He and K*He formed in the 
lowest SpH and dpH states, respectively, all vibrational 
levels including the lowest ones were found to be pop¬ 
ulated to varying extents depending on the laser wave¬ 
length [26]. Therefore we tend to discard predissociation 
of Rb*He exciplexes as the origin of the peculiar shape of 
the Rb+He transients, although we cannot strictly rule 
it out. More insight into the Rb+He dynamics may be 
provided by further measurements using electron and ion 
imaging detection. 


V. SUMMARY 

This experimental fs pump-probe study discloses the 
competing dynamics of desolvation and solvation of ex¬ 
cited and ionized states, respectively, of Rb atoms which 
are initially located at the surface of He nanodroplets. 
The generic feature of the pump-probe transients - the 
time-delayed appearance of photoions - is shown to re¬ 
sult from the falling back of ions into the droplets when 
the ionization occurs at an early stage of the desorption 
process. This interpretation is backed by the experimen¬ 
tal observation of the opposing trend when measuring 
the yield of unfragmented He nanodroplets containing a 
Rb+ ion. Furthermore, mixed quantum-classical model 


calculations based on one-dimensional pseudo-diatomic 
potentials confirm this picture qualitatively. The limited 
quantitative agreement with the experimental results is 
attributed to the use of model potentials in the calcu¬ 
lations, which do not account for the transient response 
of the He density upon excitation and ionization of the 
Rb dopant atom. Much better agreement may be ex¬ 
pected from three-dimensional time-dependent density 
functional simulations [16, 21] of the full pump-probe se¬ 
quence which are currently in preparation. 

Pump-probe dynamics similar to the Rb+ case is ob¬ 
served when detecting Rb+He molecular ions which pri¬ 
marily result from photoionization of Rb*He exciplexes. 
The peculiar structure of the Rb+He transients as well as 
extended model calculations indicate that photoassocia- 
tive ionization is an additional mechanism of forming He- 
containing ionic complexes in fs experiments. However, 
the dynamics resulting from the additional photoasso- 
ciative ionization channel cannot unambiguously be dis¬ 
tinguished from predissociation of Rb*He exciplexes in 
high-lying vibrational levels of the bpS-state. 

These results shed new light on the interpretation of 
the Rb+He pump-probe transients measured previously 
by ionizing via the lowest SpH excited state [1, 24]. The 
signal increase at short delays was interpreted as the 
manifestation of the formation dynamics of the Rb*He 
exciplex by a tunnelling process. Possibly the compet¬ 
ing desorption of the excited neutral and the fall-back of 
the photoion actually more crucially determines the rise 
time of the Rb+He signal in those transients. This issue 
will be elucidated in future experiments using two-color 
pump-probe ionization via the 5p droplet states and at 
low laser repetition rate so as to exclude concurrent ef¬ 
fects by subsequent laser pulses. Furthermore, we will 
investigate the photodynamics of metal atom-doped He 
nanodroplets in more detail by applying refined detec¬ 
tion schemes such as ion and electron imaging [14, 16] 
and coincidence detection [43, 44]. 
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